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1 Introduction 

The Belgium Priest G. Lemaitre, in 1933, was the first to study spherically 
symmetric distributions of matter without pressure in a non-static, inhomogc- 
neous cosmological model [Tj. For such a model, being dust sourced, one has 
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the Lemaitre-Tolman-Bondi metric (for spherical coordinates (r, 0, </>)) 

(dR\ 2 

ds 2 = dt 2 - ^-9zl—dr 2 - R(t, rf (d9 2 + sin 2 Odcf) 2 ) . (1.1) 



Here E{r) is an arbitrary function (a local curvature function) and R(t, r) is 
both temporal and spatially dependent. For R — the dynamical equation 
(or a first integral of the Einstein equations) 



R 2 = 2E(r) + + ±R 2 (1.2) 

R o 

holds, where M(r) is an arbitrary (mass) function and where A is a cosmological 
constant. Equation (1.2), which exhibits an energy interpretation of E(r) as 
well, was solved parametrically by Lemaitre [1] in terms of the Weierstrass 
elliptic phi-function, sigma function, and zeta function. Special case solutions 
were also obtained by R. Tolman, H. Bondi, B. Datta, as referenced in the paper 
of G. Omer [2J, for example. 

More generally the present paper, which is a continuation of work initiated 
in [3] , concerns a parametric solution of the nonlinear differential equation 

m 2 = « (i.3) 

V\t) 

with applications, where f(x) = aoa; 4 + 4aix 3 + 6a,2X 2 + Aa^x + is a quartic 
polynomial and n > is a fixed whole number. Even the case n = 1 is of some 
special interest as it includes equation (1.2), for example. Namely, for n = 1 we 
write (1.3) as 

y(tf - By(t) 2 + Ey(t) - K + A + (1.4) 

which is (1.2) for E = 0, D = 0. The case n = 2, for example, is of relevance 
in a discussion of a Bianchi V cosmological model - as we shall see. Interest- 
ingly enough, similar to Lemaitre's solution of (1.2), a solution to (1.3) is also 
expressed in terms of Weierstrass' phi, sigma, and zeta functions - as we shall 
also see. Moreover, following Lemaitre, we also express our solution in terms 
of theta functions (introduced by C. Jacobi in 1829) since the series for these 
functions are known to converge quite rapidly and are therefore amenable to 
practical calculations. 

An additional impetus for the present work is the nice paper of G. Krani- 
otis and S. Whitehouse [4] which presents the most general exact solution of 
inhomogeneous relativistic cosmology, based on the Szekercs-Szafron family of 
metrics [5j[6] (also see [7]) 

ds 2 = dt 2 - e 2 «*.W) ( dx 2 + dy 2 ) _ e 2«(t, X ,y,r) dr 2^ ^ 

where again pressureless matter is assumed, and the condition ^ ^ is im- 
posed. The other case ^ = provides for a family that generalizes Friedmann 
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and Kantowski-Sachs models [5j. The functions a, (3, by way of the Einstein 
field equations, assume the form 

e f}(t,x, v ,r) = $(t,r)e v( - r ' x > y \ 

(1.6) 

e a(t,w) = ^(r)$(t,r)^(t,x,y,r), 

where i/(r, x, y) has the form 

e -v(r,x,») = ^(^(3.2 + y 2j + 2B(r)x + 2C(r)y + D(r) (1.7) 

and where <!>(£, r) satishes 

93 9$ 9^ 

${t,r) — (t,x,y,r) = —(t,r) + $(t,r) — (r,x,y), 

4 [A(r)D(r) - B(r) 2 - C(r) 2 ] = /i(r)" 2 + A'(r), (1.8) 
d$\ 2 , \ , s 2M(r) A , N9 

Thus the third equation in (1.8), which (mathematically) is equation (1.2) is 
(again) equation (1.4) with E = 0, D = 0. 

Barrow and J. Stcin-Schabes [S], for example. Also compare the earlier paper 
[TU] of G. Covarrubias. The 1997 text [7] by A. Krasinski can be consulted for 
a rather comprehensive analysis of inhomogencous cosmological models. This 
is supplemented by the 2011 Review Article by K. Bolejko, M.-Noelle Celerier, 
and A. Krasinski |llj - the main thrust of which is a description of various cos- 
mological observations that are explainable by the models of Lemaitre-Tolman 
and Szekeres. 

Part of our interest here is to extend portions of the Kraniotis-Whitehouse 
(K-W) paper [?]. For example: (i) We define and compute Hubble and decel- 
eration parameters for our solution of (1.4); the results are expressed in terms 
of the Weierstrass phi- function p(w) (ii) Whereas D in (1.4) is zero in the K- 
W paper, we consider the case D / 0, which amounts to a constant electric 
charge contribution to the third dynamical equation in (1.8). Here we find in 
the case of a vanishing cosmological constant a new (non-parametrized) solution 
in terms of J. Lambert's W-function. In the literature this Vy-function is also 
called the product logarithm or the Omega function. Namely, it is an inverse 
function, given a particular branch, of the function f(w) = we w (iii) We con- 
sider briefly the zeros of p(w) in the general context of equation (1.4). In [1] 
an expression is derived for the Hubble parameter in the de-Sitter phase (the 
inflationary scenario) at such a zero, where the zeros are described there by the 
1982 Eichler-Zagier formula [12] . There is, on the other hand, a quite different, 
newer formula due to Duke-Imamoglu |13j which we present (for the record) 
as it may be less familiar to Physicists (iv) We consider an automorphic (or 



Solutions of (1.8) (in both cases §7 7^ and = 0) were also found by J. 
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modular) property of our solution. Such a (deep) property is referred to (but 
not explicated) in the paper [14] of E. Abdella and L. A. Correa-Borbonet, in 
the special case of a Friedmann universe - in connection with conformal field 
theory on a torus and J. Cardy's entropy formula - the torus being constructed 
from a ratio of the two periods of the Weierstrass function piw). These authors 
also take the K-W paper as impetus for their work. 

Apart from the K-W paper and an application of the solution of (1.3) to 
Bianchi V cosmology, we find in addition (v) a family of new solutions for 
a Bianchi IX model in Lyra geometry. This extends some initial work of G. 
Bag, B. Bhui, S. Das, and F. Rahaman [TS]. (vi) Some details are provided 
that indicate why the differential equation for the second moment I 2 (t) of a 
wave function of the Gross-Pitaevskii equation is also an example of (1.4), with 
A = B = 0. Thus we can compute this moment (for a stiff fluid model), 
which is essential for the dynamic correspondence set up in [T|)J [17] between 
Fricdmann-Lemaitre- Robertson- Walker (FLRW) and/or Bianchi I cosmology 
and Bosc-Einstcin condensates - a correspondence which one may regard as 
between cosmology and condensed matter - or (very generally) between a grav- 
itational and a non-gravitational system. 

Various examples are presented, together with 3 appendices, to assist in the 
navigation between theory and application. 

As is well-known, the concept of entropy in the arena of quantum statis- 
tical mechanics was developed by John von Neumann. One of his co-authors, 
and a staff member of his electronic computer project at the Institute of Ad- 
vanced Studies at Princeton during the year 1952-1953 was A. Edward (Eddie) 
Nussbaum, to whom this paper is dedicated. Professor Nussbaum was one of 
the most influential teachers and a mentor to the second named author - who 
expresses herewith his manifold appreciation for that influence and inspiration 
over the course of many years. 

2 The Solution of (1.3) in Terms of Theta Func- 
tions 

Associated with the polynomial f(x) in (1.3) are the invariants g 2 ,g3 and 
the discriminant A given by 



a Oi a 2 
Oi a 2 a 3 , 
a 2 a 3 a 4 

(2.1) 

Let ei, e2, e 3 denote the roots of 4a; 3 — g 2 x — g 3 = and fix a root x$ of f(x) = 
0, which we assume is non-repeated. If the leading coefficient ao of f(x) is 
indeed non-zero, then xq is non-repeated A ^ 0, as shown in chapter 3 of 



d&J ■ A O 2 

g 2 = a a 4 — 4aia 3 + 3a 2 , 

def- n 3 2 2 

.93 = aga 2 a4 + 2aiQ,2a 3 — a 2 — aoa 3 — 0^(24 



A = .92 - 27g 



3' 
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[TB] . for example. In [3] we have derived the following parametric solution of 
equation (1.3), where we refer to Appendix A following section 4 of this paper 
for a review of the definitions of the classical Weierstrass functions p(w) — 
p(w;uj 1 ,aj 2 ),cr(w), and ((w): 



y 



x 



/'(so) 



A[p(w + c)- /"(zo)/24]' 
f'(xo) 



X 



4[p(w + c)- f"(x )/24] 



dw + 5, 



(2.2) 
(2.3) 



for integration constants c, 5. Moreover, for a choice of any wq with p(u>o) 
/"(a;o)/24, p(wo) ^ e\, e 2 , e3 we have for n = 1,2 respectively in (2.3) 



/'(so) 
4p'(wo) 



log 



<r(tt) + c — wo) 
cr(w 4 



w 



2(w + c)((w ) 



5, 



Xof'{x ) , /'(x ) 2 p"K) 



+ 



2p'( Wo ) 16p'K>) 3 



log 



cr(u; 



w J 



/'(so) 2 
16p'(w ) ; 

+(w + c) 



[C(w + c + w ) 



p'(wo) 



C(w ) 



f C(w + 

/'(^o) 2 
16 



it('u; 
- w )] 
2p(w ) 



p'(wo) 2 



Wq) 



2p"(u» )C(w ) 
p'(w ) 3 



(2.4) 

(2.5) 
+ S, 



by formulas 1037.06, 1037.11, respectively, in [19J . In case n = 0, one has the 
(non-parametrized) Biermann- Weierstrass solution HU H2] 



!/(*) = 2/(0) 



/(tf(0))VV(t) + (p(*) 



/"fa(Q)) 

24 



f(y(Q))f""(y(Q)) 

24 



P(*) 



/"fa(o)) 

24 



f(y(Q))f""(y(o)) 

48 



(2.6) 

of (1.3). Here the periods u>i,(d2 of p(t) are constructed from the invariants 
52,93 of /(#) in (2.1), according to Appendix A. Formula (2.6) simplifies quite 
a bit in case y(0) is a root of f(x), which is the case if, for example, the boundary 
condition y(0) = is imposed: 



y(t) = y(o) 



/'(y(o)) 



p(*) 



24 



(2.7) 



which is reminiscent of equation (2.2). One can replace y(0) in (2.7) by any 
non-repeated root xq of f(x) in fact. We specifically choose the order of the 

def 

roots ej, j = 1,2,3, by setting e,- = p(uij/2;uj 1,1x2), where W3 =' u>i +u>2- 

We proceed now to a theta function expression of the solution - the impor- 
tance of which was indicated in the introduction. There is, unfortunately, no 
uniformity in the definitions/notation of theta functions as indeed considerable 
variations exist across various texts. Here we shall employ, specifically, the def- 
initions and notation of the Chandrasekharan text |18j . together with certain 



5 



results therein: For a complex number z € C and r S the upper i-plane II + 
(i.e. Irar > 0), and q d = e mT 

def. 



6{z,t) 
9i(z,t) 



def. 



def. 



def. 



2 ^](-l) n g (" + 2) 2 sin((2n + 1)ttz), 



2^] (? (»+3) 2 cos((2n + l)^), 

n=0 



l + 2^(-l)"q" 2 cos(27rnz), 

n=l 

oo 

1 + 2^q" 2 cos(27rnz). 



(2i 



For r fixed, these are entire functions of z. Obviously 0{z, r) is an odd function 
of z and the 0j(z, t), j = 1, 2, 3, are even functions of z. 

For the choice r = £ II + , the theta function expression of p(w) = 

p(w;wi,W2) is given by 



p{w) =e + — 
wf 



9'(0,r) 



,(0,r) 



e j ( 












8 ( — , 













(2.9) 



for j = 1,2,3, where 9'(z,t) denotes partial differentiation with respect to z. 
Similarly there are theta function expressions of the sigma and zeta functions 
given by 



a(w) 



CM 



UJl r niw 2 /uii. 



w 



(2.10) 



a(w) u>i u>i 





( w 




1 ° 






UJl g 


i w 











def def 

for r =' u%luj\, 7]i =' C (^i/2) • In particular 



A I U — Wq \ 

O-(u-W ) = c - AmWou/oJl \ ui ' J 
<t(u + Wo) Q ( U+WQ T \ 

V Wl ' J 



(2.11) 
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which with (2.10) gives 



a(w + c — wo) , , , 
log _)_ \ : _ °( + 2(w + c)C(wo) 



cr(w + c + w ) 



loE 



l(w + c)e' (^,r) 



(2.12) 



where in (2.12) the constant 771 no longer appears. Strictly speaking, in the 
argument for (2.12) one has log(ziZ2) = logzi + log 2:2 + 2niN, \oge z = z, for 
N G {—1,0, 1} , — 7r < Im(z) < it, for the principal branch of the logarithm. 
However we absorb the 2mN into the integration constant 5 in (2.3). That is, 
in summary we can express the solution of equation (1.4) parametrically by y 
and t in (2.2) and (2.4), or equivalently (using (2.9), (2.12)) by 



x 



24 



?'(0,r) 



9i(0,r) 



w-\-c 



far) 



XqW 



f'(xo) 
Apf (wo) 



Icjc 



, T 



2(w + c )6'[^,T 



W + C+Wq 



far) 



,(2.13) 



+ <5, (2.14) 



for j = 1,2, 3, r = U2/LU1, xq — a non-repeated root of /(x) = 0, and for the 
theta functions defined in (2.8); the ej, wq arc defined following (2.1) and (2.3). 
Here for n = 1 we have /(x) = Bx A + Ex 3 — Kx 2 + Ax + D, in accordance with 
the notation in (1.4). 

Equation (2.13) for y (or equivalently equation (2.2)) holds independently 
of whatever value n assumes. Only the expression for t (as a function of w) 
depends on n, by (2.3). For n = 2, for example, t is already given by formula 
(2.5), which we also express in terms of theta functions as follows: 



XqW + 



-Xq/'(xq) f(x ) 2 p"(w ) 

[ 2p'(w ) + I6p'(w ) 3 



ft ( w+c-\ 
. W , x , V "1 

4r?i — (w + c) + log 



+ W 



/'(xo) 5 



16wip'(wo) g 



Arji(w 



I I W + C+Wq 



W + C+Wq 



^'l 



w 'x /'(x ) . > /'(x ) 2 
+ (w + c) — — — -({wo) - 



p'(w ) 



16 



W-\-C — Wq 
UJl 



I ( W-j-C— Wq 



W-\-C — Wq 



2p(w ) | 2p"(w )((w ) 
p'(w ) 2 p'(wo) 3 



, T 



(2.15) 



by (2.10) and (2.11), where (again) 771 d =' C(wi/ 2 ) 
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Although the cases n = 0,1,2 treated here are sufficient for the physical 
applications we have in mind, one can in principle carry out the integration in 
(2.3) for n > 3 by using the recursion formula 1037.12 on page 312 of |19j . 

Naturally associated with the solution (2.2), (2.4) of equation (1.4) are the 
Hubble and deceleration parameters 

TTU ,def.y(t) def. y{t)y{t) 

H{t) = — , q(t) = — — , (2.16) 

y(t) y(t) 

respectively. Using that (2.4) was derived by integration of (2.3) - i.e. 

dw Xo+ 4[p(w + c)- f"(x )/24] [ ' 

for n = 1, one applies parametric differentiation to deduce the following formu- 
mA -4f'(x )p'(w + c) 



(4x [p(w + c) - £^2i] + /'(x )) 



«(*) 



(2.18) 



p"(w + c) \p{w + C )-^ 



4x 



p'(w + c) 



- 2 



p{w 



f"(x ) 



f'(xo) 



Here f{x ) = 4Bx 3 + 3Ex 2 - 2Kx + A, f"(x )/24 = f x% + f x - g. These 
formulas simplify in case xq = is a root of f(x) = - say (for example) D = 
in (1.4), the case considered in [4]: 



H(t) 



-4p'(w + c) 
~1W) 



(2.19) 



_ p"(w + c){p(w + c)-f"(0)/24} 
q[) ~ p'(w + c) 2 



compare formulas (48) and (50) in [J. Since p(z) satisfies the differential equa- 
tion 

p'(z) 2 =4p(z) 3 -g 2 p(z)-g 3 , (2.20) 

one can use that p"(w + c) = 6p(w + c) 2 — g 2 /2 in formulas (2.18), (2.19). Note 
that a slight misprint occurs in formula (46) of [4] for R(t). Namely, a minus 
sign is needed. Similarly formula (65) should read p'(wo) = ~i\/gl}- 

We have assumed so far that the root Xo was non-repeated. For repeated 
roots, solutions of (1.3) in fact are generally easier to obtain and they usually 
involve elementary, non-elliptic functions. As an example, take f(x) = — 4x 4 + 

2x 2 + x/y/2 + 1/16 which has xq d =' —y/2/4 as a repeated root: (x — xq) 2 is 



a factor of f(x). Equation (1.3) (with n = 0) is u'(x) 2 = -4u{x) 4 + 2u{x) 2 + 
u(x)/y/2 + 1/16, which has elementary, non-elliptic solutions 

u(x) = u T (x; 5) = — ■= ' -, (2. 21) 

V ; TV ' 2V2[ v / 2±sin(x-(S)] V ' 

as discussed in 0[23], for example. One can "deform" the solution u+(x;0) to 
obtain a solution u(x,t) of the modified Novikov-Veselov equation 

u t = u xxx + 2Au 2 u x . (2.22) 

Namely, one takes 

«(*,*) = u + (x + 2t;0) = _ _ . . (2.23) 

Note that a solution of the differential equation 

Z/i(t) 2 = Vi (t) 2 [Br + E iyi (t) m - K lVl (t) 2m + A lVl (t) 3m + D m (t) im ] , 

(2.24) 

where m ^ is a fixed whole number, can be obtained from a solution y(t) of 
equation (1.4). For this we set y(t) = yi(t)~ m . Then for B = m 2 Bi, E = m 2 Ei, 
K = m 2 Ki, A — m 2 Ai, and D = m 2 D 1 , equation (2.24) is transformed to 
equation (1.4). 

3 A Modular and Elliptic Property of the Solu- 
tion of (1.3) 

The phi- function is homogeneous of degree —2 in the sense that for any 
A € C — {0}, p(Xw; Xuji, X0J2) = X~ 2 p(w;lji,l}2)- In particular for A = w-f 1 and 
t = G II + (as in (2.10)) p{w\Wx, 0J2) = 2 p 1; rj , which means that 

the phi-function is determined by the special phi-function p(w\\;r) that we 
focus on here, and which we denote by p(w; r) for an arbitrary t € II + . p(w; r) 
satisfies the wonderful modular property [24] 



w ar + b 



d 



(cT + d) 2 p(w;T) (3.1) 



for 7 = 



a b 

c d 



€ SX(2,Z); i.e. a,b,c,d€ Z and det"/ = 1. The elliptic (or 
quasi- periodicity) property 

p{w + rar + n; t) = p(w; r) (3.2) 

is also satisfied for m, n € Z. The modular property is fundamental for the dis- 
cussion in |14| regarding solutions of the Friedmann equation (in various settings 
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such as that of a 4-dimensional radiation dominated universe, for example) and 
connections to conformal field theory (CFT) on a torus (defined by the modular 
parameter r), and with some linkage to the Cardy entropy formula. We offer a 
bit more on this later. Here A(l,r) ^ (see Appendix A) and r depends on 
the cosmological constant. 

Wc illustrate how the modular property translates to a concrete modular 
property of the solution of equation (1.3). Thus we incorporate the dependence 
of the solution on r and write y(w; r) for y in (2.2), where the choice of lattice 
in Appendix A is .if = Jz? (1, r): 

y(w; t) = x + —— { Wo A] ■ (3.3) 

4[p(w;r) - /"(x )/24] 

Here we have chosen the integration constant c in (2.2) to be zero. Then for 
a b 

7 



G SL (2, Z) with 7 • r = 7 ' %g 



77777TT- (3-4) 



The idea now is simple: In (3.4), apply (3.1) and then use (3.3) to express 
p(w; t) in terms of y(w; r). The result is 



w ar + b 



CT + d CT 
Xq 



(3.5) 

f'(x )[y(w,T) - Xq] 



i^pl [(ct + d) 2 - 1] [y(w; t) - xo] + f'(x )(cr + df ' 

In addition to the modular property (3.5) of the solution (3.3), equation (3.2) 
provides for the (less profound) elliptic property 

y(w + niT + n; r) = y(w; r) (3-6) 

for to, n G Z. 

If (again) = is a root, our formulas simplify. Suppose in fact, for 
example, that also f"(0) = which means that K = in (1.4). Then equation 
(3.5) reduces to 

( C T + d)- 2 y(w,T), (3.7) 



CT + d CT + d 

which also follows directly from (3.1), (3.3). The case K = is discussed in 
section 5 of [3], where a flat universe with non-zero cosmological constant is 
considered. One can't help but note the comparison of (3.7) with (3.1). 

The authors in |14j . seeking new connections between the Friedmann equa- 
tion and Cardy formula, set up the "chain of connections" (using their phase 
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on page 5): Friedmann equation — > Weierstrass equation — > torus — >■ CFT 
partition function Z(t) — > Cardy formula, where for q d =' e 27rlT 

Z(t) = trace q L °~ c / 24 = trace e~ m (3.8) 

with c = a central charge, H = Lq — c/24 a Hamiltonian corresponding to the 
Virasoro generator Lq, (3 = —2-kit an inverse temperature, and with the Cardy 
entropy formula [25) given by 

Here Lq also denotes the eigenvalue of Lq and, for simplicity, the anti-holomorphic 
sector (which involves r) is not considered. As mentioned above r depends on 
the cosmological constant A; it also depends on the spacetime dimension. 

We are positioned to mimic this chain of connections as follows. The Fried- 
mann equation is replaced, more generally, by equation (1.3). The Weierstrass 
equation (or solution) is given by (3.3) (and by (2.3)). r in (3.3) defines a cor- 

def. 

responding complex torus T =' C/Jz?(1,t), which (as is well known, based on 
(2.20)) is an elliptic curve parametrized by the phi-function and its derivative. 
The remainder of the chain is as before, given by (3.8), (3.9). 

John Cardy's great insight was that the key property of modular invariance 
of CFT partition functions was deeply consequential. In particular it lead to 
his asymptotic density of states formula and hence to his entropy formula 
(3.9) being a special case of the general formula that also incorporates the non- 
holomorphic sector. 

If we write (3.1), (3.2) as 

p(— ^-jJ^-T^) = e 2 "^(cT + d) k p(w,T), 
\CT + d CT + d J 

(3.10) 

p{w + mr + n;T) = e - 2n]{m2r+2mw) p{w;T) 

for j, k £ Z (namely for j = 0, k = 2), we see that p(w; r) is a (mcromorphic) 
Jacobi form of index j = and weight k — 2 [24] . Similarly, given (3.6), (3.7), 
we can regard y(w; r) as a Jacobi form of index and weight —2. Strictly 
speaking, a growth condition at infinity is also imposed in the definition of 
a Jacobi form (which we shall not have a concern with here) - a condition 
analogous to that imposed in the definition of a modular form [26j . The elliptic 
genus of an N = (2, 2) Super CFT with central charge c = 6j is, for example, a 
Jacobi form of index j and weight k = 0. 
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4 Some Applications 



Wc present some examples to provide further clarity and applications of the 
formulas in section 2. 

Example 1. Among various modifications of Einstein's GR was one intro- 
duced by G. Lyra in 1951 in which a gauge function was introduced to affect 
Riemannian geometry that rendered it a generalization of Weyl's 1918 geom- 
etry. More specifically, the purpose of the gauge function was to remove the 
non-integrability of vector length under parallel transport. Many cosmological 
and string models based on Lyra geometry have been studied. In particular the 
Bianchi IX model in this geometry was studied in [15]. Here, for the metric 

ds 2 = -dt 2 + a(tj 2 dx 2 + b(t) 2 dy 2 + [b{t) 2 sin 2 y + a(t) 2 cos 2 y] dz 2 

-2a(t) 2 cosydxdz (4.1) 

and a constant flat potential V((f>) = 2A, the modified Einstein equations yield 
the relation <f> = ^jab 2 for an integration constant <j)Q, and the assumption 
a = b" leads moreover to the first integral 

b 2 = * - * + D lb - 2 - 2 , (4.2) 

— 1 2(n z — n) 

for n 7^ 0,±1 and an integration constant D\. The authors in [15| found solu- 
tions of (4.2) only for £>i = 0, n = 2, |, |, and |. A solution for n — 2, D 1 ^ 
(say D\ = To) was found in [3]. We indicate here how to construct an infinite 
family of solutions with A^O corresponding to 

, , def 3A 4 + 4A 3 , 4 1 ± V2i , s 

I>i=i3i(A) ^ - 2 ,A^0,-1,--,— y~, (4.3) 

where for A = 1 we obtain the solution just mentioned, corresponding to D\ = 

7_ 

12- 

For n = 2, a = b 2 and equation (4.2) is transformed to the equation 

d(t) 2 = -a(t) 2 + L(t) + ^, (4.4) 

which is equation (1.4) for B = -1, E = f, K = 0, A = 0, D = 4Di. From 
(2.1), .92 = -4Di, .93 = -4Di/9 = fla/9, and A = -16D 2 (41?! + |) . A = 
=> D\ = 0, — jL, which for Z?i = -Di(A) in (4.3) means that A = 0, — | and 
3A 4 +4A 3 = -1. The latter equation is = 3A 4 +4A 3 + 1 = (A + 1) 2 (3A 2 -2A + 1), 
which gives A = -1, (1 ± y/2i)/3. That is, by (4.3), A ^ for D 1 = L»i(A) so 
by the remarks following (2.1) f(x) = —x A + |x 3 +4Di(A) = has no repeated 
roots, where we note that xq = —A is indeed a root. As /'(—A) = 4A 2 (A + 1) 
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and /"(— A)/24 = — A(3A + 2)/6, we obtain the parametric solutions 

A 2 (A + 1) 



a = -A + 



p(w + c) + A(3A + 2)/6' 

(4.5) 



t = -Aw + 



X 2 (\ + l) 



cr(w + c-w ) „, x>/ x 
cr(w + c + u> ) 



(5 




of equation (4.4), by formulas (2.2), (2.4), for A ^ 0, -1, -f , (1 ± y/2i)/3. Then 
in the metric (4.1), b(t) 2 = a(t). 

Example 2. Consider next, briefly, the metric 

ds 2 = -dt 2 + X(t) 2 dx 2 + e 2bx [Y(t) 2 dy 2 + Z(tfdz 2 \ (4.6) 
for a Bianchi V anisotropic cosmological model, where 6 ^ is a constant. 

R(t) [_X"(t)Y(£)i?(i)] 1/ ' 3 , as a consequence of the Einstein equations, satisfies 
a differential equation 

R(t) 2 = f(R(t))/R(t) 4 (4.7) 

for a suitable 4 th degree polynomial f(x) with leading coefficient b 2 . A key, 
useful observation here is that the quantity /function 

de d- (4.8) 

X(t)Y(t) X(t)Z(t) Y{t)Z{t) 
' X(t)Y(t) ~ X(t)Z(t) ~ Y(t)Z(t) 

is actually a constant, independent of t, which we may therefore denote by Dq. 
Dq is used, in fact, to construct the constant term of f(x). Other coefficients 
of f(x) involve radiation and matter constants, and the gravitational constant, 
where we assume that the energy momentum tensor is given by that of a perfect 
fluid. For full details of these remarks and further related information the reader 
can consult [H [27] . 

Since equation (4.7) assumes the form (1.3) with n = 2, R(t) is given para- 
metrically by y,t in equations (2.2), (2.5) - or cquivalently by (2.13), (2.15). 
On the other hand, the "scale factors" X(t),Y(t), Z(t) can be determined from 
R(t) by the formulas in [3j[27], and thus the metric (4.6) is also determined by 
R(t). 

Example 3. The influence of an electro-magnetic field on the collapse of dust is 
considered briefly by Krasinski in section 7 of [55] . A more extended discussion 
and list of references is found in his text [7]; see section 2.14, for example. Here, 
again with spherical symmetry assumed, the third equation in (1.8) (or, equiv- 
alently, the Lemaitre equation (1.2)) is supplemented by a term — (Q/$>(t, r)) 2 , 
where Q is a constant electric charge: 
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It will turn out (see (2.1)) that 



.92 = --7 



K{rf 
12 



53 



AQ 2 K(r) K{rf AM(r) ; 



18 



216 



12 



(4.10) 



A 



A 2 Q 6 AQ A K{r) 2 Q 2 K(r) 4 AQ 2 K(r)M(r)' 



27 18 

K{rfM(rf 3AM(r) 4 
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In particular A = when the cosmological constant A vanishes, in which case 
one expects non-elliptic solutions of (4.9). Indeed for A = 0, Shikin has 
found parametric solutions of (4.9) in terms of elementary, non-elliptic functions. 
On the other hand, however, since (4.9) is a special case of (1.4) (with B = 



E = 0,K = K(r), A = 2M(r), D 



f{x) = Ax i -K(r)x 2 +2M(r)x-Q 2 



from which one computes (4.10)) we have that its general parametric solution 
(for arbitrary A) is given already by equations (2.2), (2.4): 



4> 



x 



XqW 



A 3 _ K(r) „ 
3 ^0 2 



M(r) 



(4.11) 



A 3 _ Kjfjxo 
3^0 2 



MM 



+2{w + c)C(wq)] + S, 



log 



a(w + c — u>o) 
a(w + c + Wq) 



where now the integration constants c = c(r), 8 = 8{r) depend on r; throughout 
it is best to consider r fixed. Also, as usual, wo is a choice such that p(u>o) = 



f"(x )/24: = (2Ax 2 - K(r))/12 ^ the roots ei,e 2 ,e 3 of 4a; 3 - g 2 x - g 3 = 0. In 
general (unless Q = 0, as in [4]) the non-repeated root Xq of f(x) = in (4.11) 
is non-zero, which means (in particular) that the formulas in (2.18) for the 
Hubble and deceleration parameters do not simplify as in (2.19). Xo of course 
also depends on r. Formulas (2.13), (2.14) provide, in addition, theta function 
expressions of the solution (4.11). Although Q does not appear explicitly in 
(4.11), the solution there is indeed Q— dependent as xq also depends on Q, and 
the invariants 32,53 depend on Q as well by (4.10). 

Section 5.2.2 of [4] contains the intriguing remark that interesting physics 
and mathematics arise when equation (2.20) is solved at the zeros zo of p(z): 



P'(^o) 5 



"S3- 



(4.12) 



Appendix C contains a description of these zeros (that supplements that given in 
[4]), where (as in section 3) the choice of lattice _£f (1, r) is made for r £ H + that 
we shall also employ here. Particularly in section 5.2.2, a Euclidean universe 
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with g 2 = 0, (73 7^ 0, A < is considered. Given the presence of the electric 
charge Q, we can slightly extend part of the discourse there as follows, where 
for convenience we write K, M for K(r),M(r). Note first that for g 2 = 0, K 2 = 
4AQ 2 =s> K 3 = AAKQ 2 =s> 

AQ 2 K 4AKQ 2 AM 2 A ( M 2 2KQ 2 \ 
93 = + -M6 12" = " 3 {'A —) (4 - 13) 

by (4.10). If we take w — z$ — c in (2.18) and denote by to the corresponding 
value of t (where we could take c = if we wished) then using (4.12) we see 
that 

H(t ) = - 4 ^°H±^) (4.14) 
(-4xof"(xo)/2A + f'(x )) 2 

which by (4.13) simplifies as 
H(t ) 

for m>o 



(4.15) 

( -K 2 x 2 a Kxq i\ , , — 

for_Q#o \JKnr + — - l ) 8KQ 2 jl 

( k*o* 5Kxo ,A 2 V 9M2 V3' 

\12MQ 2 6M ' L ) 

where in the latter equation we use again that K 2 = 4AQ 2 for 172 = 0. In case 
Q = 0, the 2 nd equation here reduces to H(to) = T^J~§ since we can then choose 
xq = as in [J], where we take the speed of light = 1. In fact by (4.10) since 
also K = 0, A = _3A!^1 < for A ^ (i.e. A f 0), which shows that x = 
is a non-repeated root. In general we see that for A > 0, H (to) is a multiple of 
the Hubble constant y/A/3 for the de Sitter model with inflationary scale factor 




Aexp (yf*) 



for some constant A. 
Kraniotis and Whitchouse discuss quite a bit more (again for the case Q = 
0) regarding other possibilities also such as A < (which leads to periodic 
solutions), A > with 172 7^ 0, (73 = (for a non-Euclidean universe), or 
(generally) A 7^ 0, g% 7^ 0, 33 7^ - including some "bouncing" models, and 
asymptotic inflationary models. 

Equation (4.9) with Q — also appears in the work of N. Meures and M. 
Bruni on ACDM cosmology [30], for example. The general solution (4.11) of 
course is not given there. In fact these authors also take the curvature constant 
K = 0, and thus they solve the equation 

• o 2M A 9 . 

$ 2 = + -$ 2 . (4.16) 

$ 3 y J 
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In general, note that the equation 

m 2 - 

has the non-elliptic solution 



A 

W) 



By{t? 



1/3 r 



sinh ( -^VBt 



-i 2/3 



(4.17) 



(4.18) 



for example. 

Consider now the special case A = 0. We present solutions of (4.9) not found 
in Shikin [29]. In this case f(x) = —K{r)x 2 + 2M(r)x — Q 2 is quadratic. Hence 

the roots of f(x) — are given by xq = (m ± M 2 — KQ 2 ^j /K, where for con- 
venience we write M, K for M(r), K(r); we also assume that K ^ 0. There are 2 
possibilities to think about: (i) M 2 ^ KQ 2 , (ii) M 2 = KQ 2 . If M 2 ^ KQ 2 , the 

roots are distinct (non-repeated). Thus for xq 

equations —Kxq 
for the solutions 



M ± \JM 2 - KQ 2 ) /K the 
M = tV^ 2 ~ KQ 2 , A = 0, with equation (4.11) provide 



<I> 



M ± \J M 2 - KQ 2 yjM 2 - KQ 



K 



T 



2 [p(w 



K/12] 



(4.19) 



M ± JM 2 - KQ 2 JM 2 - KQ 2 
W =F 



K 

-2(w + c)((w )}+5 



2p'(w ) 



log 



a(w + c — w ) 
<j(w + c + wq) 



for p(w ) = -K/12 ^ ei, e 2 , e 3 . 

To obtain quite different solutions we therefore focus on the 2 nd possibility: 
M 2 = KQ 2 , in which case xo = M/K is a repeated root: f(x) = —K{x — 
M/K) 2 , and the general formulas (2.2), (2.4) that lead to (4.11) do not apply. 
However, equation (4.9) now assumes the form 



<j> 2 



-K($ - M/K) 2 /<5> 2 , 



(4.20) 



which has (non-parametric) solutions in terms of Lambert's W- function men- 
tioned in the introduction: 



(4.21) 



where B ^ is an integration constant that depends on r. Note that since 
M{r) 2 = K{r)Q 2 and K(r) / Owe also have M(r) ^ for Q ^ 0, and 
Q 2 /M(r) = M(r)/K(r). One can derive (4.21) from the general fact that the 
differential equation 

x{t) = b+^ (4.22) 



a 
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with a, b ^ has the solution 



x(t) 



-b 



W(Be 



-tb 2 /a\ 



(4.23) 



for B ^ since W{z) satisfies W'(z) = W{z)[z{\ + W(z))}- 1 . 

Apart from the application here, the Lambert VF-function has been applied 
in a variety of disciplines ranging from statistical mechanics and quantum chem- 
istry to enzyme kinetics, the engineering of thin films, and the physiology of 
vision 1311. 



Example 4. A dynamic correspondence between FLRW and/or Bianchi I cos- 
mology and Bose-Einstein condensates (BECs) governed by a time-dependent, 
harmonic trapping potential was set up in |17j . A cosmological constant A^ 
was present, where the spacetime dimension d > 3 was arbitrary. Thus an 
extension of work of James Lidsey [16] was realized. The correspondence is 
presented in Tables I and II below where cosmological parameters (scale fac- 
tors, scalar pressure and energy density p^,, Hubble parameters) are matched 
with wavepacket parameters expressed in terms of the harmonic trapping fre- 
quency uj(t) and moments Ij(t), j = 2,3,4 (with Izit) > 0) of a wave func- 
tion of the Gross-Pitaevskii (G-P) equation. Here t is "laboratory" time that 
one obtains from the "cosmic" time in the Einstein field equations. We also 
assume constancy of the atomic interaction parameter in the G-P equation. 
In Table I, a(t) is the scale factor for the FLRW cosmological model and 

H(t) d = a(t)/a(t) is the corresponding Hubble parameter. In Table II, R(t) = 
(Xi(t)X 2 (t) ■ ■ ■ Xd-iit)) 1 /^ 1 *) is the average scale factor with Xj(t) the scale 

factor in the j th spatial direction and Hn{t) d =' R{t)/R(t). Also we have 

def 

Kd = 8irGd for Gd the gravitational constant. 



Table 1: BEC O FLRW correspondence 



h ^ 
11/ Ah 

[(d-l)(d-2)I 4 -A d ]/K d o 

[{d - 2)uj 2 I 2 - \d- l)(d - 2)h + Ad] /K d <-> 



a 

2{aH) 
H 2 
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Table 2: BEC <->■ Bianchi I correspondence 



h 
h 
It/Ah 



(d-2) 



A, 



(d-2) 2 j (d-2) j 



(d-iy 



(d-1) 

A 



/K d 
+ K. 



/K d 



fl 2(d-l) 

2(rf-l)(i?( d - 1 )i/ fl ) 

(d-l) 2 #l 



-2(d 



P4> 

l)K d D/{d 



2) 



The Tables are based on an equation of state = (7 — l)p<j, with 7 > 0. 
7 = 6/(d — 1), for example, corresponds to a stiff perfect fluid. The moments 
satisfy the conservation law 

2I 2 (t)h(t)-I$(t)/4= a constant d = X. (4.24) 

The two basic equations derived in |17j in conjunction with Table I (with the 
help of (4.24)) were 

2 _ ao7(d-l) I\ _ 2a 2A d I 2 

^ z [7(d-l)+2]/2' 4 / [ 7 (d-l)-2]/2 + 2 ) 1 J 

for an integration constant ao- These equations govern the time-dependent 
trapping frequency u)(t), and hence they also govern the external potential 
V(r,t) = uj(t) 2 r 2 /2. In particular, for a stiff perfect fluid the second equation 
here is written 

for a = 8ao, which is equation (1.4) with A = B = there. Instead of the para- 
metric solution given by (2.2), (2.4), we have an alternate (simpler) parametric 
solution 

h = ap(w;g 2 ,g 3 ) +b 

(4.27) 

t = -a 2 ((w;g 2 ,g 3 ) + abw + 5, 

as constructed in [3], where a 3 = (d - l)(d - 2)/2A d , b = X(d- l)(d - 2)/6A d , 
,g 2 = 2a(d - l)(d - 2)A 2 /3A d , ff3 = 2A 3 (d - l) 2 (d - 2) 2 /27A 2 - a. 

The Bicrmann-Weierstrass solution (2.6) applies directly to provide for an 
explicit solution of the scale factor a{rj) of one of the two the Friedmann equa- 
tions, where eta is "conformal" time - the other equation being a local conser- 
vation of energy equation. This is discussed in [32j [33], for example. There 
energy and matter in the universe are assumed to be a perfect fluid consisting 
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of radiation, non-relativistic matter, and a cosmological constant. The modular 
property (3.1) and the "chain of connections" of [14] discussed in section 3, and 
the solution (4.11) in the special Friedmann case (with Q = 0, xq — 0) are 
also discussed in |37| . for example, where a holographic description of the early 
universe is considered. 



Appendices 

A Definition of the Weierstrass phi, sigma and 
zeta functions 

Given the central importance of the Weierstrass phi function p(w) for the 
present work we recall briefly, for the reader's convenience, its construction/ 
definition. A more detailed account is available in [T51 |2"T1 153] . 

Let oji,u>2 be non-zero complex numbers. Since the imaginary parts of 
a non-zero complex number z and its reciprocal are related by Im z~ l = 
— (Im z) \z\ 2 , one has that Im oj 2 /ui\ ^ if and only if Im lji/lo 2 ^ 0. In 
particular we assume that Im ui 2 /u>i > 0, which is equivalent to the assumption 
Im loi/u>2 < 0. The corresponding lattice Jz? = j£?(wi, co 2 ) generated by u)i and 
u>2 is defined to be the set of points ui = mu)\ + nu>2 where m and n vary over 
the set of whole numbers. The lattice Jzf gives rise to the phi function 



p(w) = -j+ J2 
w L — ' 

ugif-{0} 



1 1 

(w — Ld) 2 UJ 2 



(A.l) 



which is also denoted by p(w;Jf), or by p(w;u>i 1 u!2). p{w) is a meromorphic 
function, which is doubly periodic with periods uj\, ui2- Thus, by definition, p(w) 
is an elliptic function. p(w) has double poles at w — uj S , and it satisfies the 
differential equation 

p'(wf = 4p(w) 3 - g 2 (ui,uj2)p(uj) - g 3 (oj 1 ,u 2 ) (A.2) 

for invariants 

52(^,^)^60 ^ g 3 (0Ji,cJ 2 ) d ^U0 ]T ± (A.3) 



where, moreover, 



A(u>x,u; 2 ) d ^ g 2 (uJi,uJ 2 ) 3 - 27g 3 (u 1 , u> 2 ) 2 + 0. (A.4) 
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Conversely, it is an amazing fact that if two complex numbers g 2 and 53 
are given that satisfy the condition g\ — TIg\ ^ 0, then there exists a pair of 
non-zero complex numbers wi,W2 with Im W2/W1 > such that 92(^1, UI2) = fj2 
and 53(^1,^2) = 53i for g 2 ( w ii ^2) and 53 (wi, u> 2 ) defined in (A. 3) with respect 
to the lattice Jz? = Jz?(wi,u; 2 ) generated by uj 1 and w 2 . Thus from g 2 and 33 
one can also construct the corresponding phi function p(w; 102) (according 
to definition (A.l)), which in this case we also denote by p(w; g 2 , g 3 )- 

Associated with p(w) are the Weierstrass sigma and zeta functions a(w) and 
((w), respectively: 

C'(w) de J- - P ( W ), lim (CM - -) de d- 0, 

(A.5) 

P^H <ie/. , . g(w) de/. 

; ; = CW, hm = 1. 

a{w) w^o w 

B Theta function notation 

The notation for Jacobi theta functions varies widely and wildly from text to 
text, as we indicated in section 2. There, in (2.8), we used the notation (and def- 
initions) presented in chapter 5 of K. Chandrasekharan [T5]. On the other hand, 
many researchers (including Kraniotis and Whitehouse |3]) employ the notation 
in the Abramowitz-Stegun Handbook [M]. Thus for the reader's convenience 
we set up the following comparison: 

K. Chandrasekharan Abramowitz-Stegun 
6{z,t) Bifaq) 
0i(z,t) d 2 (irz,q) 
6 2 (z,t) 9 A (irz,q) 
6 3 {z,t) 9 3 {nz,q) 

for q = e 7 ™ 1 ", r € Also the notation in Whittaker- Watson [21] is the same 
as that in Abramowitz-Stegun, where the 84 (z, q) in [21] is initially denoted by 
6(z,q). 

The periods wi, w 2 of the phi-function in Appendix A are denoted by 2w, 2a/ 
in [34]: lui = 2w, uj 2 = 2uj'. Applications of theta functions to FLRW cosmology 
also appear in [36], where the notation of [21] is employed. 



C On the zeros of p(w; r) 

The physical significance of the zeros of p(w; r) has been pointed to in Example 
3. For the sake of completeness we provide a description of these zeros. 

With definitions (A. 3), (A. 4) in mind, we introduce the normalized Eisen- 
stein series 

W = /, / i S 2 (1,t), E 6 (r) de J- 1^.93(1, r), (C.l) 
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and normalized discriminant 

def. A(l,r) 



A„(r) 



27 A(l,r) 



64 2 tt 12 64 x 1728 tt 12 



(C.2) 



for t £ II + . Then the zeros Zo of p(u>;r) are given by the following explicit 
integral formula of M. Eichler and D. Zagier [T2l |2"4"] : 



zq = m + — + jit ± 



27TI J T E 6 (<j) 3 / 2 



(C.3) 



for m, n <= Z, where the integral is over the vertical line er = xo + (yo + t)i, t > 0, 
in n + commencing at t = xo + j/o*j J/o > 0. 

An alternate formula for the zeros was found by W. Duke and O. Imamoglu 
[13] . who "dcuniformized" the Eichler-Zagier formula and expressed the zeros 
in terms of the classical modular invariant 



j(T)^-E4rr/A n (r). 



(C.4) 



The generalized hypergeometric series 



p Fq(ai,...,a p ;bi,...,b q ;z) =' ^ 



de/. ^ (0l)„ • • • (0 p )„ « T 



(£>l)n ' ' ' (6g)n «! 

are needed, where |z| < 1 and each (bk) n 7^ for the Pochhammer symbol 
(o)„ d = J ' a(a + l)(a + 2) ■ ■ ■ (a + n - 1) = F( f + , n > 1, 



(C.5) 



r(o) 



(C.6) 



(a) = 1. 



Forp = 2, q = 1, 2-Fi(ai, «2; &i; z) is the standard Gauss hypergeometric function 
F(ai,a 2 ; h; z). Let 



de/. -i-v/6 , . , 1728 

C2 = —Z , U = U{r) = l-— r-r- 

3tt j(rj 



1728A„(r) 

£ 4 (r) 3 : 



(C.7) 



and choose the principal branch of it 1 / 4 . Then the zeros of p{w; r) are ±zo 
where zq is given by the Duke-Imamoglu formula 



1+T 



2(1 



C 2 U 



1/4 _3 



^2(5, 3, l; 4, |; u) 



2 fi(^,^,l;l-u) 



(C.8) 



u is given in terms of r in (C.7). r, conversely, is given in terms of a as a 
quotient of hypergeometric functions: 



r (A) r (ts) 2 ^i (A' il; !; 1 _ 



-i + 



(C.9) 
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In a completely different context, formula (C.8) also has an application (in- 
terestingly enough) in the work of Conte, Grundland, and Huard [35], for exam- 
ple, on isentropic ideal compressible fluid flow. More specifically, (C.8) is used 
to conclude the boundcdness (and hence the physical relevance) of a certain 
elliptic solution of theirs. 

Note that by (C.l), (C.2), (C.4), one also has the following expression for j: 

i M 1 79sMkL>! de J- 1728ff 2 (l,r) 3 (cm) 

j(t) = 1728—— — - = — - — -3 . — rj, (L.10) 

A(l,r) 2 (1, rf - 2753(1, ry 

which is the form expressed in [3] , where a flat universe with cosmological con- 
stant A ^ and j(r) = is considered (as we have discussed in Example 3), 
and where a non-Euclidean universe is also considered with j(r) = 1728. A 
full discussion of the modular invariant j(r) and of Eisenstein series and modu- 
lar forms, including explicit Fourier expansion formulas, is presented in |26j . In 
particular, we discuss there the key role that the Fourier expansion of j(r) plays 
in the duality between 3-dimensional pure gravity with a negative cosmological 
constant and extremal conformal field theories with central charge 24k, for a 
positive integer k. 
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